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Magnetic acceleration of ultra-relativistic GRB and AGN jets 
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We present numerical simulations of cold, axisymmetric, magnetically driven relativistic outflows. The outflows are 
initially sub-Alfvenic and Poynting flux-dominated, with total— to-rest-mass energy flux ratio up to ~ 620. To study 
the magnetic acceleration of jets we simulate flows conflned within a funnel with rigid wall of prescribed shape, which 
we take to be z oc r° (in cylindrical coordinates, with a ranging from 1 to 2). This allows us to eliminate the numerical 
dissipative effects induced by a free boundary with an ambient medium. We find that in all cases they converge to a 
steady state characterized by a spatially extended acceleration region. For the jet solutions the acceleration process 
is very efficient - on the outermost scale of the simulation more than half of the Poynting flux has been converted 
into kinetic energy flux, and the terminal Lorentz factor approached its maximum possible value {Too — A*). The 
acceleration is accompanied by the coUimation of magnetic fleld lines in excess of that dictated by the funnel shape. 
The numerical solutions are generally consistent with the semi-analytic self-similar jets solutions and the spatially 
extended acceleration observed in some astrophysical relativistic jets. In agreement with previous studies we also find 
that the acceleration is significantly less effective for wind solutions suggesting that pulsar winds may remain Poynting 
dominated when they reach the termination shock. 

Keywords: jets acceleration; gamma-ray bursts; active galactic nuclei; pulsar wind nebulae. 



1. Introduction 

There is strong evidence for relativistic motions in 
jets that emanate from active galactic nuclei (AGNs). 
The Lorentz factors of blazar jets lie in the range 
~ 5-40 [1-3]. 

In the case of AGNs there have indeed been in- 
dications from a growing body of data that the as- 
sociated relativistic jets undergo the bulk of their 
acceleration on scales that are of the order of those 
probed by very-long-baseline radio interferometry. In 
particular, the absence of bulk-Comptonization spec- 
tral signatures in blazars has been used to infer that 
jet Lorentz factors ^ 10 are only attained on scales 
^ 10" cm [4]. 

The first theoretical clues to the necessity of rela- 
tivistic motion in GRB's arose from the compactness 
problem [5] . The requirements that the source be op- 
tically thin can be used to obtain direct limits on the 
minimal Lorentz factor, F > 100 [6, 7]. Recently, su- 
perluminal expansion was observed in the afterglow 
emission of GRB 030329 [8]. 

The main source of power of AGN and GRB jets 
is the rotational energy of the central black hole [9, 
10] and/or its accretion disk. The naturally occur- 
ring low mass density and hence high magnetization 
of black-hole magnetospheres suggests that the rela- 



tivistic jets originate directly from the black-hole er- 
gosphere as Poynting-dominated outflows, whereas 
the disk surface launches a slower, possibly non- 
relativistic wind that surrounds and confines the 
highly relativistic flow. Although, dissipative pro- 
cesses may directly transfer the electromagnetic en- 
ergy to emitting particles, the commonly held view 
is that it is first converted into the bulk kinetic en- 
ergy and only subsequently channeled into radiation 
through shocks and other dissipative waves [11-16]. 
In the limit of ideal MHD such conversion can be 
achieved only via magnetic forces and the efficiency 
of this mechanism is the main subject of our investi- 
gation. 

Since most of the acceleration takes place far 
away from the source the space-time is basically flat. 
We also use an isentropic equation of state p = Qp^ , 
where Q =const and s = 4/3. Since we are inter- 
ested in the magnetic acceleration of cold flows, we 
make Q very small, so the gas pressure is never a dy- 
namical factor. This relation enables us to exclude 
the energy equation from the integrated system and 
overcome the stiffness of relativistic MHD in mag- 
netically dominated regime. 

Given the condition of axisymmetry the poloidal 
magnetic field is fully described by the so-called mag- 
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netic flux function ^{r,z), the total magnetic flux 

enclosed by the axisymmetric loop circle r, z =const. 
Moreover, for stationary flows there are 5 quantities 
that propagate unchanged along the magnetic fleld 
lines and thus are functions of ^ alone. These are 
k, the rest- mass energy flux per unit magnetic flux; 
li, the angular velocity of magnetic field lines; /, the 
total angular momentum flux per unit rest-mass en- 
ergy flux; /i, the total energy flux per unit rest-mass 
energy flux; and Q, the entropy per particle. For cold 
flows {Q = 0, w = pc^) we have k = pTvp/Bp, 
Clr = - VpB^/Bp, I = -I/2iTkc + ru^, and 
/i = r (1 + (t) , where F is the Lorentz factor, Vp is 
the magnitude of the poloidal velocity, Bp and B'^ 
are the magnitudes of poloidal and azimuthal mag- 
netic field, / = crB'^/2 is the total electric current 
flowing through a loop of cylindrical radius r, a is 
the ratio of the Poynting flux to the matter (kinetic 
plus rest-mass) energy flux, and Per = —Vtl /2'Kk(? is 
the Poynting flux per unit rest-mass energy flux. It is 
easy to see that the Lorentz factor F cannot exceed 
11. 
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Fig. 1. Evolution of a along the magnetic field line ^I" = 
Q.S^max in models with a = 620 (solid line), 310 (dashed 
line), 155 (dash-dotted line), 78 (dotted line) and 39 (dash- 
triple-dotted line). 

2. Numerical Setup 

In the simulation we deal with winds, or unconfined 
flows, and jets, or flows confined with a funnel. In 



order to avoid complications due to numerical diffu- 
sion through the jet boundary we use solid funnels 
of paraboloidal shape, z oc r° , where z and r are the 
cylindrical coordinates. 
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Fig. 2. Fct (upper branch) and F (lower branch) along the 
magnetic field line ^ = Q.S^max (solid lines), along the mag- 
netic field line "I* = 0.5'I'm,ax (dashed lines), and along the 
magnetic field line 4' = Q.2^rnax (dash-dotted lines) in model 
(TO = 620. 

The initial configuration corresponds to a non- 
rotating, purely poloidal magnetic field with approx- 
imately constant magnetic pressure across the fun- 
nel. The plasma density within the funnel is set to a 
small value so that the outflow generated at the inlet 
boundary can easily sweep it away. In order to speed 
this process up the longitudinal component of veloc- 
ity inside the funnel and at the inlet is set to 0.7 c. 
In fact, we use the same type of initial and boundary 
conditions as described in our study of lower mag- 
netization jets [16]. We constructed a grid of models 
with different funnel power geometry (a = 1,3/2,2) 
including unconfined wind, different initial magneti- 
zation (cr = 10^600), and both solid and differential 
rotation at the base. 

3. Results 

Figure 1 shows the evolution of a along the jet 
boundary for the models with a = 3/2, solid rota- 
tion, and the maximum magnetization at the base 
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ao = 39, 78, 155, 310 and 620 (for solid rotation ao in- 
creases with r). One can see that aU sohitions exhibit 
transition to the particle dominated regime {a < 1). 
Prior to reaching the equipartition they are described 
by the same law a oc {r/ric)~^^^ and after this the 
evolution of a slows down significantly. The depen- 
dence of 'the equipartition' radius req on ctq can be 
approximated by 

req ^ 0.079(7^ V(c if 10 < (70 < 620. (1) 

Figure 2 confirms similar evolution inside the jets. 
The normal pressure at the jet boundary is close to 

Pn cc R~^, where R is the spherical radius. In astro- 
physical conditions this pressure has to be balanced 
by the pressure of confining medium. 

Figure 3 shows our results for the unconfincd 
wind. One can see that magnetic acceleration in this 
case is much less effective in agreement with previous 
studies [17, 18]. 

4. Application to GRB and AGN Jets 

and PWN 

The initial energy-to-mass flux ratio of jets in our 
simulations yields an upper limit on the terminal 
Lorentz factor T^o = A*. In order to make further 
comparisons of our numerical models with observa- 
tions we need to select suitable dimensional scales. 
The key scale in the problem of magnetic accelera- 
tion is the light cylinder (or the Alfven surface) ra- 
dius, ric- If the jets are launched by a rapidly rotating 
black hole in the center of an GRB or AGN then 

ric - 4rg = 2 X 10*^(^/73 Mq) cm, 

where Vg = GM/c^. In this estimate we assume that 
the angular velocity of the magnetic field lines is half 
of that of a maximally rotating (rotation parameter 
a ~ 1.0) black hole. According to the results shown 
in equation (1), the jets enter the matter-dominated 
regime at a cylindrical radius 

Teq ~ 5.9 X lO^(At/15O)2-i(M/3M0) cm. 

The corresponding distance from the black hole is 

Req ~ 5.9 X 10i°(Ai/150)2-i(M/3Mo)(0.1/ej) cm, 

where Qj is the jet opening half-angle. Strong shock 
waves can appear only if cr < 1. For GRB jets, with 
fjL > 150 and M ~ 3Mq this distance is remark- 
ably close to the radius of Wolf-Raye stars. Thus, 



in the coUapsar scenario the conditions for develop- 
ment of strong fast shocks and associated gamma-ray 
emission are already satisfied when the jets break 
away from the progenitor star. The actual location 
of shocks, however, also depends on the scale of the 
central engine variability which puts it to larger dis- 
tances. 

Lower Lorentz factor of AGN jets imply lower 

magnetization parameter, /i < 30. Combined with 
the black hole mass M ~ IO^Mq this gives Req — 
2 X 10^'^cm which is similar to the size of the "blazar 
zone" inferred from the observations. Thus, the slow 
nature of magnetic acceleration of relativistic jets is 
in very good agreement with the observational con- 
straints. 

The velocity profile of jets with a solid body ro- 
tation at their base shows strong tangential discon- 
tinuity at the jet boundary which favors the photon- 
breading mechanism of gamma-ray emission from 
AGN jets [19-21]. On the contrary, the gradual de- 
cline of Lorentz factor towards the jet boundary 
which we see in the models with differential rota- 
tion at the base significantly reduces its efficiency. 
This implies that gamma-gay observations can be 
used to determine whether the AGN jets originate 
directly from the black hole magnetospheres or from 
the magnetospheres of their accretion discs. 

Like other researches [17, 18] we find that mag- 
netic acceleration of highly magnetized unconfined 
winds is less ineffective. This is a rather uncom- 
fortable conclusion since the successful MHD model 
of Pulsar Wind Nebulae [22-24] requires the pulsar 
wind to be particle dominated at the location of its 
termination shock. The two possible solutions to this 
problem suggested so far utilize the alternating struc- 
ture of magnetic field in the wind from oblique ro- 
tators ('striped wind', [25]). Firstly, the alternating 
magnetic field can dissipate before reaching the ter- 
mination shock, e.g. [26] , with eventual convertion of 
the released heat into the bulk motion energy. Sec- 
ondly, it can dissipate inside the shock layer of the 
termination shock itself [27]. In this case one cannot 
directly use the standard shock equations in order to 
determine the downstream state. 

5. Conclusion 

In validating the basic features of the simplified semi- 
analytic solutions, our numerical results go a long 
way toward establishing an "MHD acceleration and 
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Fig. 3. Unconfined wind solution. Left panel: Lorentz factor along three different magnetic field lines: solid line - * = O.S^max, 
dashed line - = O.bipmax, dash-dotted line - * = Right panel: Ta (solid line), n (dashed line) and T (dash-dotted 

line) along the magnetic field line with 5' = O.S^max as a function of cylindrical radius. 



coUimation paradigm" for relativistic astrophysical 
jets. In particular, they demonstrate that even jets 
with extremely high initial magnetization can be ef- 
fectively accelerated via the ideal magnetic mech- 
anism with more than a half of the Poynting flux 
converted into the bulk motion energy of the flows. 
The highest Lorentz factor reached in the simulations 
is r = 300 which is well within the range deduced 
for GRB jets. The slow character of magnetic accel- 
eration allows dissipationless energy transport over 
large distances, the property which is deduced from 
observations and which is rather difficult to explain 
in other models of jet generation. 

Interested reader can find more details about the 
method of our simulations and results for AGN jets 
in [16]. 



References 

[1] Cohen M. H. et al, ApJ, 658, 232, (2007). 

[2] Jorstad S. G. et al, AJ, 130, 1418, (2005). 

[3] Jorstad S. G. et al., ApJS, 134, 181, (2001) 

[4] Sikora M. et al., ApJ, 625, 72, (2005). 



[5] Rederman M., Ann. N.Y. Acad. Set, 262, 164, 
(1975). 

[6] Krolik J.H. et al., ApJ, 373, 277, (1991). 

[7] Piran T., Phys. Rep., 314, 575, (1999). 

[8] Taylor G.B. et al., ApJ, 609, LI, (2004). 

[9] Blandford R. D.et al, MNRAS, 179, 433, (1977). 

[10] Lovelace R. V. E., Nature, 262, 649, (1976). 

[11] Begelman M. C. et al.. Reviews of Modem Physics, 

56, 255, (1984). 

[12] Blandford R. D. et al., MNRAS , 169, 395, (1974). 

[13] Vlahakis N. et al, ApJ, 596, 1080, (2003a). 

[14] Vlahakis N. et al, ApJ, 596, 1104, (2003b). 

[15] Beskin V.S. et al., MNRAS, 367, 275, (2006). 

[16] Komissarov S. S. et al., MNRAS, 380, 51, (2007). 

[17] Beskin V.S. et al, MNRAS, 299, 341, (1998). 

[18] Bogovalov S. et al, MNRAS, 305, 211, (1999). 

[19] Aharonian F. et al. Science, 314, 1424, (2006). 

[20] Derishev E.V. et al. Physical Review D, 68, 043003, 
(2003). 

[21] Stern B.E. et al., MNRAS.tmp.lWiS, (2007). 

[22] Kenel C.F. et al, ApJ, 283, 694, (1984). 

[23] Komissarov S. S. et al, MNRAS, 349, 779, (2004). 

[24] Del Zana L. et al, A&A, 421, 1063, (2004). 

[25] Michel F. C., Reviews of Modem Physics, 54, 1, 
(1982). 

[26] Arons J., astro-ph/0710. 1326B, (2007). 

[27] Petri J. et al, A&A, 473, 683, (2007). 



